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ABSTRACT 

Texture  analysis  is  performed,  on  multibeam 
sonar  signal  returns  discriminated  angularly 
by  beamforming.  A  collection  of  fourteen  tex¬ 
ture  features  are  computed  via  co-occurrence 
matrices  and  data  reduction  is  then  performed 
using  a  principal  components  transformation. 
Acoustical  boundaries  (boundaries  between 
regions  with  homogeneous  acoustioal  proper¬ 
ties)  are  evident  from  the  features.  Results 
indicate  that  seafloor  bottom  characteristics 
can  be  extracted  from  these  texture  features. 

SEAFLOOR  IMAGERY  FROM 
MULTIBEAM  SONAR  DATA 


nnc  quality  inspected  9 


The  impetus  for  generating  seafloor  acoustic 
imagery  from  multi  be?™  sonar  datnis  that . 
more  detailed  information  is  available  from  the 
returning  pulse  than  is  discernible  through  sim¬ 
ple  computation  of  a  set  of  bottom  depths  (de 
Moustier,  1986).  This  is  due  to  the  additional 
information  contained  in  the  signal  intensity  and 
to  the  typically  high  hydrophone  sampling  rates 
which  provide  finer  across-track  resolution  than 
is  possible  through  beamforming.  The  resulting 
imagery  displays  the  response  of  the  seafloor  to 
the  acoustic  pulse,  which  is  a  combination  of  the 
local  topography  as  well  as  the  bottom  composi¬ 
tion.  Past  efforts  have  attempted  to  provide 
methods  for  segmenting  acoustic  imagery  into 
geoacoustic  provinces  via  various  image  texture 
processing  methods  to  allow  automated  classifi¬ 
cation  of  the  seafloor  (Reed  IV  and  Hussong, 
1989;  Bourgeois  and  Walker,  1991).  However, 
the  process  of  generating  imagery  from  the  raw 
beamfomied  data  filters  at  a  low-pass  rate  and 
distorts  the  image,  and  full  resolution  acoustic 
imagery  is  too  noisy  for  successful  fine-scale 
segmentation.  In  this  work  analysis  is  performed 
on  the  returning  energy,  as  a  function  of  time, 
that  has  been  angularly  discriminated  by  beam¬ 
forming.  This  approach  allows  analysis  on  data 
that  have  not  suffered  the  additional  generation 
of  noise  and  introduction  of  smoothing  that 
occurs  when  forming  an  image. 

The  data  used  are  from  the  Navy’s 
Sonar  Array  Survey  System  (SASS),  and  com¬ 
plete  navigational  information  is  not  available; 
we  do  not  have  information  as  to  the  actual  ori¬ 
entation  of  the  sonar  array  from  ping  to  ping 
Since  the  data  are  not  fully  georeferenced,  each 
ping's  swath  is  processed  individually.  The 
results  of  the  texture  feature  extraction  for  each 


ping  is  then  appended,  line  by  line,  to  form  a 
multiband  texture  feature  set.  Parameters  used 
for  georeferencing  can  be  used  later  to  find 
accurate  correspondences  between  seafloor 
location  and  bottom  characteristics  as  comput¬ 
ed  from  the  texture  analysis. 

MULTIBEAM  DATA  PROCESSING 

System  Hardware 

Multibeam  bathymetric  sonar  systems 
such  as  SASS,  the  SeaBeam  series,  and  the 
EM  100  are  capable  of  collecting  data  which, 
after  proper  processing,  may  be  used  to  accu¬ 
rately  map  the  bottom  of  the  ocean.  It  is  possible 
to  obtain  both  bathymetry  and  sidescan-like 
images  fF©m  the  original  raw  data.  Information 
about  texture  of  the  ocean's  bottom  may  also  be 
generated  directly  from  this  data. 

For  the  SASS,  the  sonar  energy  from 
the  projector  array  located  on  the  hull  of  the  sur¬ 
vey  ship  impinges  on  the  bottom  of  the  ocean  as 
a  narrow  beam  perpendicular  to  the  ship’s  head¬ 
ing.  The  echo  from  this  swath  is  received  by  an 
array  of  hydrophones  mounted  athwartships 
(perpendicular  to  the  projector)  under  the  sur¬ 
vey  vessel.  The  beam  width  depends  on  the 
angle  of  arrival  to  the  bottom  and  increases  with 
inerr  -sing  angle.  The  unsteered  beam  width  is 
0.7";  at  45°  it  is  0.9°. 

Data  Processing 

Processing  starts  by  reading  the  raw- 
data  recorded  by  SASS  during  a  survey  mission. 
A  detailed  description  of  the  original  data  is 
given  by  Bourgeois  (1991)  and  one  of  the  com¬ 
plete  processing  is  given  by  Kaminsky  et  al. 
(1992).  The  beamforming  process  is  done  next 
by  performing  delay,  filtering,  and  summing 
operations  through  the  use  of  Fourier  trans¬ 
forms.  This  yields  an  array  of  return  intensi¬ 
ties — an  intensity  for  each  beamformer  bin  and 
each  sample  time— allowing  the  returning  ener¬ 
gy  to  be  resolved  into  angular  bins.  Each  of  the 
K  beamformer  bins  k  correspond  to  a  particular 
steering  angle  <bk. 

The  beamfomied  data  then  give  us  a 
time  history  of  the  energy  received  from  the  K 
look  directions.  This  data  must  be  further  pro¬ 
cessed  to  determine  the  time  corresponding  to 
the  center  of  the  beam,  tr  The  peak  of  the  enve¬ 
lope  corresponds  to  the  intersection  of  the  maxi¬ 
mum  response  axis  (MRA)  of  the  steered  beam 
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and  the  area  ensonified.  Once  the  time  tc  has 
been  determined  by  one  of  many  methods  avail¬ 
able  (Kaminsky  et  al.,  1993)  the  bathymetry  of 
the  area  surveyed  can  be  obtained.  That  is,  exact 
points  ( x,y,z )  that  determine  the  source’s  posi¬ 
tion  are  obtained  from  each  (<}>*,  tj  pair. 

The  intensity  data  produced  by  the  beam- 
former  may  also  be  used  to  produce  sidescan¬ 
like  images  that  show  the  backscatter  intensity. 
One  method  of  forming  image  pixels  is  by  sum¬ 
ming  the  total  energy  returned  at  a  given  sample 
time.  In  this  way,  one  port  and  one  starboard 
pixel  are  generated  for  each  sampling  instant. 
Due  to  the  beamforming  process,  the  image 
quality  can  be  improved  by  discarding  the 
returning  energy  in  directions  not  associated 
with  the  dominant  pulse  for  the  given  sample 
time  using  a  spatialAemporal  window. 

In  this  paper  the  energy  returned  along 
each  beam  is  examined  as  a  function  of  time. 
Specifically,  the  amplitude  variations  in  the 
returning  pulse  for  each  steered  beam  are  ana¬ 
lyzed.  These  amplitude  variations  are  descriptive 
of  the  roughness  of  the  seafloor  in  the  region 
illuminated  by  the  acoustic  pulse.  If  this  region 
is  small  enough  so  that  the  sediment  properties 
may  be  considered  homogeneous  then  the 
amplitude  variations,  or  signal  ‘texture,’  are 
descriptive  primarily  of  the  physical  roughness 
of  the  seafloor. 

Variation  in  local  slope  is  a  large 
source  of  ambiguity  in  acoustic  imagery:  drasti¬ 
cally  different  backscatter  intensities  are 
obtained  depending  upon  the  acoustic  incident 
angle,  making  automated  analysis  of  this  form  of 
data  extremely  difficult.  The  rest  of  this  paper  is 
devoted  to  texture  feature  extraction,  manipula¬ 
tion,  and  its  interpretations. 

TEXTURE  ANALYSIS 

exture  is  a  property  of  pixel  values  and  their 
spatial  relationships.  It  is  insufficient  to 
describe  pixel  values  without  considering  their 
spatial  relationships  to  other  pixel  values.  To 
quantitatively  analyze  textures,  we  must  use 
measures  that  are  not  only  a  function  of  individ¬ 
ual  pixel  values  but  also  multiple  pixel  values 
and  their  relative  spatial  positioning. 

There  has  been  a  lot  of  work  done  on 
the  analysis  of  textures.  The  methods  usually  fall 
into  certain  categories,  for  example,  co-occur¬ 
rence  statistics  (Haralick  ct  al.,  1973);  power 
spectrum  methods  (Lendaris  and  Stanley,  1970); 
use  of  Markov  models  (Cross  and  Jain,  19S3; 
Manjunalh  and  Chellappa,  1991);  structural  anal¬ 
ysis  (Haralick,  1979);  and  fractal  analysis 
(Pentland,  19S-1). 

Conners  and  Harlow  compare  different 
texture  algorithms  and  conclude  that  the  co¬ 
occurrence  methods  (which  they  refer  to  as  the 


Spatial  Gray  Level  Dependence  Method)  is  'he 
most  powerful  method  because  of  its  ability  to 
discriminate  among  a  set  of  textures  (Conners 
and  Harlow,  1980).  Mastin  et  al.  used  co-occur¬ 
rence  methods  for  SAR  (synthetic  aperture 
radar)  imagery  of  water  (Mastin  et  al.,  1985). 
Consequently,  we  decided  to  use  an  analysis 
based  on  co-occurrence  statistics  because  of  the 
performance  of  the  method  and  also  the  nature 
of  the  data  acquisition  (proper  georeferencing 
information  was  not  available).  The  data  avail¬ 
able  to  the  authors  do  not  allow  georeferencing 
or  even  referencing  between  pings.  Therefore, 
we  only  analyze  the  texture  along  the  one¬ 
dimensional  path  of  each  ping. 

Co-Occurrence  Statistics 

A  histogram  is  an  estimate  of  the  first 
order  statistics  of  an  image  (or  a  region).  The 
normalized  histogram  is  computed  as 

P(t)  =  ^,i=0,  1, . . .  2*  -  1  (1) 

—where  N(i)  is  the  number  of  pixels  in  the  image 
(region)  with  intensity  value  i,  N  is  the  total 
number  of  pixels  in  the  image  (region),  and  b  is 
the  number  of  bits  per  pixel  in  the  image. 

The  analog  to  the  histogram  for  second 
order  statistics  is  the  co-occurrence  matrix.  The 
co-occurrence  matrix  is  also  computed  in  a 
“census”  fashion  by  counting  pairs  of  occur¬ 
rences  of  pixels  values  given  a  certain  spatial 
relationship  for  the  pair.  The  normalized  co¬ 
occurrence  statistics  are  computed  as 

P(t„  i2,  d,  6 )-  -  '2I  =  Did,  6 )  (2) 

for  pairs  of  pixels  at  locations  x,  and  x2  having 
values  i,  and  i2>  respectively.  The  distance  mea¬ 
sure  D(d,0)  states  that  the  spatial  relationship  of 
the  pair  of  pixels  is  that  they  are  'ocated  at  a  dis¬ 
tance  magnitude  d  apart  and  at  an  angle  0  (or 
9+tt)  from  each  other. 

A  complete  set  of  co-occurrence  statis¬ 
tics  would  cover  all  values  of  d  and  0  over  a 
meaningful  range.  The  values  for  0  would  vary 
between  0  and  it  using  some  number  of  discrete 
steps.  The  values  for  d  would  range  from  1  up  to 
some  distance  where  the  correlation  between 
pixels  is  still  significant. 

In  practice,  several  co-occurrence 
matrices  a. e  computed  for  several  pairs  of  (d,0). 
Figure  1  shows  several  computed  co-occurrence 
matrices  for  a  simple  example  2  bit/pixel  image. 

One  of  the  disadvantages  of  the  method 
of  co-occurrence  matrices  is  the  potentially  large 
amount  of  data  computed  for  different  pairs  of 
d  and  0.  Only  four  of  the  many  possible  co¬ 
occurrence  matrices  are  computed  in  Figure  1. 
However,  co-occurrence  statistics  are  powerful 
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FIGURE  1.  Sample  co-occurrence  calculations  (Haralick  and  Shapiro,  1992). 
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in  that  they  are  invariant  to  monotonic  intensity 
transformations  (Haralick  and  Shapiro,  1992). 

Texture  Features 

Texture  features  can  then  be  computed 
from  the  co-occurrence  matrices.  Ballard  and 
Brown  (1982)  suggest  a  set  of  five  features,  but  a 
more  comprehensive  set  of  14  features  is  given 
by  Haralick  et  al.  (1973).  It  is  often  not  clear  how 
these  features  relate  to  observable  phenomena. 
However,  they  have  been  demonstrated  to  be 
useful  in  classifying  images.  Using  the  co-occur¬ 
rence  statistics,  these  14  features  are  computed 
for  each  beamformer  bin.  In  this  application, 
only  a  co-occurrence  matrix  for  d=l  and  0=0  is 
computed.  We  keep  0=0  since  the  pings  are  not 
georeferenced.  The  distance  d  is  kept  small 
since  we  expect  the  primitives  to  be  relatively 
small.  However,  we  can  still  compute  more  sets 
of  these  features  for  various  values  of  d,  but  d=  1 
for  initial  results. 

We  will  follow  the  notation  used  in 
Haralick  et  al.  (1973): 

•  p(i,j)  is  the  (i,j)'h  entry  in  the  normalized 
co-occurrence  matrix, 

•  px(i)  is  the  Ith  entry  in  tire  marginal-probabili¬ 
ty  matrix  obtained  by  summing  the  rows  of 
co-occurrence  matrix. 


•  Ng  is  the  number  of  gray  levels  in  the  image, 

•  Py  U)  =  2  i  P  O'.J). 

•  Pr. y  (*)  =  2  f=l  2  %  1  P  (’-T)-  k  =  = 

2, 3,  -,  2 Nr 

•  Pt-yiV  =  2  ii\  2  Nj^  t 

*  =  lt-j!=0,  -  1. 

•  e  is  a  small  constant  to  prevent  taking  the 
log  of  zero. 

The  features  computed  are  as  follows: 

1.  angular  second  moment:  a  measure  of 
homogeneity  of  the  image, 

/i  =i  (3) 

i=i  i 

2.  contrast:  a  measure  of  the  contrast  or 
the  amount  of  local  variations  present, 

-V  £  \ 

fi  =  Y^\X^pU.f)l  (4) 

inti  i=i  ,=i 

3.  correlation:  a  measure  of  intensity  linear 
dependencies  in  the  image, 


f  .S/J-iSjvi  (uMU) 

It  I T 


where  /ij.  p j,  a,  and  ny  arc  the  means 
and  standard  deviations  of  p ,  and  pv. 
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4.  sum  of  squares  variance:  a  measure  of 
the  variation  in  the  image, 

M2)P(i-  3)-  (6 

5.  inverse  difference  moment: 

N9  %  1 

-4=22  i  +  (7 

<=i  v 

6.  sum  average: 

2N, 

fa  =  X  *?*♦»(*)■  (8 

i=2 

7.  sum  variance: 


/?  =  £(* 
i-2 


These  features  are  computed  for  each 
ping  along  all  256  bins.  Actual  processing  leaves 
off  approximately  30  bins  on  each  end  since 
these  represent  “near  horizontal"  directions  that 
produce  no  useful  return  information  and  are 
mostly  noise.  Also  bins  with  “dropouts"  have 
been  discarded.  We  know  that  some  pings  are 
missing  and  that  the  ones  we  have  are  not  neces¬ 
sarily  parallel  to  each  other.  Due  to  lack  of  posi¬ 
tioning  data,  we  display  pings  (rows  in  the 
image)  as  if  they  were  parallel  to  one  another. 
Thus,  the  features  computed  for  each  ping  are 
appended  to  one  another  forming  14  “texture 
feature  images"  that  form  the  multidimensional 
feature  set. 


8.  sum  entropy: 


DATA  REDUCTION 


2N, 

far  “2  Px*y(i)log  +  €].  (10) 

i= 2 

9.  entropy: 

/9=  -2  2p(u  j)  togWO)  +  «!•  (ii) 

_ i  l  j- 1 

10.  difference  variance: 

fio  =  variance  of  px_y.  (12) 

11.  difference  entropy: 

/ll=  -  2  Ar-yC)  lOglAr-yO  +  «)-  (13) 

i-Q 

12.  information  measure  of  correlation  1: 
for  this  and  the  next  feature,  HX  and 
HY  are  the  entropies  of  px  and  py, 
respectively  and 

HXY=  -2  2  P (W)  log(p(t,y)+  e), 

»=1 /*f 

//AY1  =  -2  itpifJ)  log|px(i)Pj,0)  +  «)>  and 

i=i  >i 


It  is  difficult  to  interpret  large  data  sets. 

Therefore,  if  we  can  reduce  the  amount  of  out¬ 
put  to  the  essentially  “useful  part"  of  the  data, 
the  reduction  will  make  it  easier  to  manipulate 
the  data  and  interpret  it. 

Data  reduction  is  also  performed  on 
the  mumucaiu  Lj.ni i_v  u u_: li ic  data  itself 
(Kaminsky  et  al.,  1992).  However,  in  this  case  we 
are  concerned  with  the  reduction  cf  the  output 
from  the  texture  analysis. 

We  have  a  fourteen-dimensional  fea¬ 
ture  set  that  was  produced  from  the  texture 
analysis.  We  cannot  just  indiscriminately  throw 
away  some  of  the  texture  features  since  we  do 
not  know  which  ones  will  be  useful  in  the  char¬ 
acterization  of  the  seafloor  bottom.  In  order  to 
pull  out  only  the  most  useful  features,  we  need 
to  find  the  principal  components  of  the  feature 
space.  This  will  give  us  a  feature  space  where 
each  of  the  new  “features"  are  decorrelaled. 

Also,  the  importance  of  each  new  feature  basi¬ 
cally  corresponds  to  the  relative  size  of  its  asso¬ 
ciated  eigenvalue.  So,  we  now  also  have  a  way 
of  distinguishing  those  components  with  the 
greatest  information  content. 


HXY2=  -2  2  PxWPyO)  logMOPyO')  +  el. 


.  i  j--i 


(  HYX-HYXX 
J[2  ~  max | //A',  HY]' 


(14) 


13.  information  measure  of  correlation  2: 
fa  =  (1  -  exp(-2.0(//AY2  -  //AT)))'*.  (15) 

14.  maximum  correlation  coefficient: 

/,,  =  (Second  largest  eigenvalue  of  Q)l/2,  (10) 


where, 


k  P,(0l>fk) 


Principal  Components 

The  feature  space  is  currently  fourteen- 
dimensional.  We  can  reduce  the  dimensionality 
of  the  feature  space  significantly  by  first  trans¬ 
forming  the  data  into  a  new  feature  space  (still 
of  the  same  dimensionality).  This  new  feature 
space  should  be  one  in  which  the  data  between 
the  features  is  uncorrelated  and  also  one  in 
which  most  of  the  “useful  information"  is  con¬ 
tained  in  just  a  few  of  those  features.  The 
Hotelling  transform  (the  discrete  formulation  of 
the  Karhunen-l/Oeve  transform)  is  used  in  this 
case  to  achieve  the  principal  components  of  the 
feature  space.  Unlike  many  other  transfonna- 
tims,  the  Hater.ir.fi  transfuim  is  data  dependent. 

The  Hotelling  transform  is  computed  as 
follows  (Gonzalez  and  Woods.  1992):  Consider  a 
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multispectral  image  consisting  of  n  bands  of  size 
N  x  N.  We  form  the  column  vector 


ond  highest  variance,  and  so  on.  We  can  recon¬ 
struct  the  original  set  of  vectors  x  using 


x  =  [XiX2 . . .  xn)'  (17) 

for  each  pixel  in  the  image.  Therefore,  each 
location  in  the  image  is  a  vector  x  consisting  of 
the  n  pixel  values  from  each  of  the  n  bands  at 
that  particular  location.  There  are  N2  such  vec¬ 
tors  in  the  multispectral  image.  The  mean  vector 
and  covariance  matrix  of  x  are  defined  as 

=  £{xj  (18) 

and 

Cx  =  FRx  -  mx)(x  -  mx);j.  (19) 

We  can  estimate  both  of  these  statistics  using 


and 

«  i 

c«=tf?2x*x'*-nvnx'.  (21) 

*=1 

Once  we  have  calculated  these  statistics  from 
the  multispectral  data,  we  compute  the.Hotelline 
transform  using  the  equation 

y  =  A(x  -  mx)  (22) 

where  A  is  the  matrix  formed  from  a  sorted  set 
of  eigenvectors  of  the  covariance  matrix  Cx 
From  this  we  see  two  important  properties: 

1.  the  transformed  data  y  has  zero  mean. 
That  is, 

m}.  =  E(y)  =  0  (23) 

2.  and  that  the  data  are  decorrelated 

Cy  =  E((y  -  m))(y  -  my)'|  =  AC^'  (24) 

which  is  a  diagonal  matrix.  The  matrix  A  is  com¬ 
posed  of  the  eigenvectors  of  Cx  and  it  will  thus 
diagonalize  it.  Therefore,  the  diagonal  elements 
of  the  covariance  matrix  of  the  transformed  data 
are  nothing  more  than  the  set  of  eigenvalues: 

\2 

Cy  =  \i  (25) 


Assuming  that  the  eigenvalues  arc  sort¬ 
ed  from  highest  to  lowest  (and  the  matrix  of 
eigenvectors  is  correspondingly  sorted),  then 
the  set  of  principal  components  as  given  in  (22) 
will  he  sorted  so  that  the  first  component  is  the 
component  with  the  highest  variance,  the  sec¬ 
ond  component  is  the  component  with  the  sec- 


x  =  A-‘y  +  mx  =  A'y  +  mx.  (26) 

We  note  that  A-1  =  A'  since  the  matrix 
is  orthonormal.  In  any  case,  suppose  we  decide 
that  the  principal  components  associated  with 
the  lowest  eigenvalues  were  of  little  use.  This  is 
often  true  since  the  magnitude  of  the  smaller 
eigenvalues  are  negligible  when  compared  with 
the  others.  Therefore,  the  images  associated 
with  these  lowest  eigenvalues  have  relatively 
small  variance,  that  is,  there  is  relatively  little 
useful  information  contained  in  these  compo¬ 
nents.  We  can  reduce  the  dimensionality  of  the 
feature  space  by  discarding  these  components. 
Suppose  we  want  to  keep  A'  of  the  n  compo¬ 
nents,  thus  reducing  our  feature  space  to  only  K- 
dimensional  instead  of  rc -dimensional.  We  can 
do  this  by  forming  a  matrix  A^  composed  of  the 
K  eigenvectors  corresponding  to  the  K  largest 
eigenvalues.  We  can  alter  (26)  to  form  approxi¬ 
mations  of  the  original  data  as 

“  ~  x  =  A'*y  +  mx7  (27) 

The  mean  square  error  between  x  and 

x  is 

(28) 

J-l  3- I 

The  Hotelling  transform  is  optimal  for 
minimizing  the  quantity  (Gonzalez  and 
Woods,  1992). 

For  the  purposes  of  feature  extraction, 
we  can  loosely  interpret  this  method  as  a  means 
to 

1.  find  a  new  set  of  features  (the  principal 
components),  such  that  these  features 
are  decorrelated  and  have  a  “measure  of 
relative  importance"  associated  with 
them,  and 

2.  reduce  the  dimensionality  of  the  feature 
space  providing  data  compression. 

RESULTS 

igures  2  through  5  show  the  first  four  princi¬ 
pal  components.  The  contrast  has  been 
improved  so  that  the  characteristics  of  the  fea¬ 
tures  are  more  visible.  The  horizontal  band 
across  the  center  of  each  feature  corresponds  to 
the  nadir  and  near-nadir  (almost  vertical)  beam 
steering  angles.  This  band  represents  a  uniform 
texture  which  is  expected  for  this  region  since, 
at  high  incident  angles,  the  acoustic  return  is 
largely  due  to  specular  reflection  instead  of 
backscatter.  The  returning  pulse  will  typically  be 
relatively  smooth  for  high  angles,  and  a  texture 
measure  of  this  sort  will  not  reveal  interesting 
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'  FIGURE  2.  First  principal  component. 


FIGURE  4.  Third  principal  component. 


FIGURE  3.  Second  principal  component. 


FIGURE  5.  Fourth  principal  component. 


features  until  the  acoustic  return  is  predomi¬ 
nately  non-spocular. 

It  is  evident  from  these  features  that 
there  are  regions  with  similar  texture  character¬ 
istics  The  texture  is  computed  from  the 


hackscatter  returns  and  represents  spatial  rela¬ 
tionships  of  the  way  in  which  the  hackscatter 
varies  along  the  pings.  Thus,  regions  with  similar 
features  represent  regions  of  similar  acoustical 
backscatlering  properties.  The  acoustical 
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boundaries  are  the  locations  of  significant  tran¬ 
sition  between  acoustical  properties. 

Figure  6  (see  the  outside  back  cover  of 
this  issue)  shows  the  first  principal  component 
(the  one  with  the  most  variance)  mapped  to  a 
terrain  obtained  from  the  bathymetric  data.  The 
center  band  running  through  the  middle  repre¬ 
sents  the  nadir  and  near-nadir  positions  as 
explained  earlier. 

CONCLUSION 

he  texture  analysis  provides  a  description  of 
the  directional  backscatter  returns  using 
second  order  statistical  properties  and  is  a  use¬ 
ful  measure.  This  is  especially  useful  even  for 
non-calibrated  systems.  Regions  of  homoge¬ 
neous  acoustical  properties  are  brought  out 
through  this  texture  analysis.  Variations  in 
backscatter  are  due  to  changes  in  bottom  types, 
surface  orientation,  and  roughness.  Due  to  the 
method  used  in  this  analysis  the  effect  of  orien¬ 
tation  is  reduced  and,  assuming  a  small  enough 
acoustic  footprint,  the  bottom  type  may  be  con¬ 
sidered  homogeneous.  Thus  these  texture  mea- — 
sures  provide  an  indication  of  micro-roughness 
(i.e.,  the  small-scale  surface  roughness  that 
exists  below  the  resolution  of  the  resolvable 
bathymetry). 

Clearly,  proper  georeferencing  infor¬ 
mation  would  improve  the  results  so  that  a  more 
accurate  classification  of  the  acoustical  proper¬ 
ties  and  the  corresponding  bottom  location 
would  result.  Thus,  the  texture  analysis  con¬ 
tributes  and  plays  an  important  role  in  seafloor 
bottom  classification  and  mapping.  In  addition, 
work  is  currently  being  performed  at  Tulane 
University  and  the  Naval  Research  Laboratory  to 
classify  the  acoustical  properties  using  neural 
models  trained  by  the  texture  features. 
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